\(\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx\) [127]
Optimal result
Integrand size = 33, antiderivative size = 189 \[
\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx=\frac {2 a^2 (39 A+34 B) \tan (c+d x)}{45 d \sqrt {a+a \sec (c+d x)}}+\frac {2 a^2 (9 A+10 B) \sec ^3(c+d x) \tan (c+d x)}{63 d \sqrt {a+a \sec (c+d x)}}-\frac {4 a (39 A+34 B) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a B \sec ^3(c+d x) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{9 d}+\frac {2 (39 A+34 B) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}
\]
[Out]
2/105*(39*A+34*B)*(a+a*sec(d*x+c))^(3/2)*tan(d*x+c)/d+2/45*a^2*(39*A+34*B)*tan(d*x+c)/d/(a+a*sec(d*x+c))^(1/2)
+2/63*a^2*(9*A+10*B)*sec(d*x+c)^3*tan(d*x+c)/d/(a+a*sec(d*x+c))^(1/2)-4/315*a*(39*A+34*B)*(a+a*sec(d*x+c))^(1/
2)*tan(d*x+c)/d+2/9*a*B*sec(d*x+c)^3*(a+a*sec(d*x+c))^(1/2)*tan(d*x+c)/d
Rubi [A] (verified)
Time = 0.56 (sec) , antiderivative size = 189, normalized size of antiderivative = 1.00, number of
steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.152, Rules used = {4103, 4101, 3885, 4086, 3877}
\[
\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx=\frac {2 a^2 (9 A+10 B) \tan (c+d x) \sec ^3(c+d x)}{63 d \sqrt {a \sec (c+d x)+a}}+\frac {2 a^2 (39 A+34 B) \tan (c+d x)}{45 d \sqrt {a \sec (c+d x)+a}}+\frac {2 (39 A+34 B) \tan (c+d x) (a \sec (c+d x)+a)^{3/2}}{105 d}-\frac {4 a (39 A+34 B) \tan (c+d x) \sqrt {a \sec (c+d x)+a}}{315 d}+\frac {2 a B \tan (c+d x) \sec ^3(c+d x) \sqrt {a \sec (c+d x)+a}}{9 d}
\]
[In]
Int[Sec[c + d*x]^3*(a + a*Sec[c + d*x])^(3/2)*(A + B*Sec[c + d*x]),x]
[Out]
(2*a^2*(39*A + 34*B)*Tan[c + d*x])/(45*d*Sqrt[a + a*Sec[c + d*x]]) + (2*a^2*(9*A + 10*B)*Sec[c + d*x]^3*Tan[c
+ d*x])/(63*d*Sqrt[a + a*Sec[c + d*x]]) - (4*a*(39*A + 34*B)*Sqrt[a + a*Sec[c + d*x]]*Tan[c + d*x])/(315*d) +
(2*a*B*Sec[c + d*x]^3*Sqrt[a + a*Sec[c + d*x]]*Tan[c + d*x])/(9*d) + (2*(39*A + 34*B)*(a + a*Sec[c + d*x])^(3/
2)*Tan[c + d*x])/(105*d)
Rule 3877
Int[csc[(e_.) + (f_.)*(x_)]*Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)], x_Symbol] :> Simp[-2*b*(Cot[e + f*x]/(
f*Sqrt[a + b*Csc[e + f*x]])), x] /; FreeQ[{a, b, e, f}, x] && EqQ[a^2 - b^2, 0]
Rule 3885
Int[csc[(e_.) + (f_.)*(x_)]^3*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_), x_Symbol] :> Simp[(-Cot[e + f*x])*(
(a + b*Csc[e + f*x])^(m + 1)/(b*f*(m + 2))), x] + Dist[1/(b*(m + 2)), Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m*
(b*(m + 1) - a*Csc[e + f*x]), x], x] /; FreeQ[{a, b, e, f, m}, x] && EqQ[a^2 - b^2, 0] && !LtQ[m, -2^(-1)]
Rule 4086
Int[csc[(e_.) + (f_.)*(x_)]*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*(B_.) + (A_))
, x_Symbol] :> Simp[(-B)*Cot[e + f*x]*((a + b*Csc[e + f*x])^m/(f*(m + 1))), x] + Dist[(a*B*m + A*b*(m + 1))/(b
*(m + 1)), Int[Csc[e + f*x]*(a + b*Csc[e + f*x])^m, x], x] /; FreeQ[{a, b, A, B, e, f, m}, x] && NeQ[A*b - a*B
, 0] && EqQ[a^2 - b^2, 0] && NeQ[a*B*m + A*b*(m + 1), 0] && !LtQ[m, -2^(-1)]
Rule 4101
Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)]*(csc[(e_.) + (f_.)*(x_)]*(
B_.) + (A_)), x_Symbol] :> Simp[-2*b*B*Cot[e + f*x]*((d*Csc[e + f*x])^n/(f*(2*n + 1)*Sqrt[a + b*Csc[e + f*x]])
), x] + Dist[(A*b*(2*n + 1) + 2*a*B*n)/(b*(2*n + 1)), Int[Sqrt[a + b*Csc[e + f*x]]*(d*Csc[e + f*x])^n, x], x]
/; FreeQ[{a, b, d, e, f, A, B, n}, x] && NeQ[A*b - a*B, 0] && EqQ[a^2 - b^2, 0] && NeQ[A*b*(2*n + 1) + 2*a*B*n
, 0] && !LtQ[n, 0]
Rule 4103
Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*
(B_.) + (A_)), x_Symbol] :> Simp[(-b)*B*Cot[e + f*x]*(a + b*Csc[e + f*x])^(m - 1)*((d*Csc[e + f*x])^n/(f*(m +
n))), x] + Dist[1/(d*(m + n)), Int[(a + b*Csc[e + f*x])^(m - 1)*(d*Csc[e + f*x])^n*Simp[a*A*d*(m + n) + B*(b*d
*n) + (A*b*d*(m + n) + a*B*d*(2*m + n - 1))*Csc[e + f*x], x], x], x] /; FreeQ[{a, b, d, e, f, A, B, n}, x] &&
NeQ[A*b - a*B, 0] && EqQ[a^2 - b^2, 0] && GtQ[m, 1/2] && !LtQ[n, -1]
Rubi steps \begin{align*}
\text {integral}& = \frac {2 a B \sec ^3(c+d x) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{9 d}+\frac {2}{9} \int \sec ^3(c+d x) \sqrt {a+a \sec (c+d x)} \left (\frac {3}{2} a (3 A+2 B)+\frac {1}{2} a (9 A+10 B) \sec (c+d x)\right ) \, dx \\ & = \frac {2 a^2 (9 A+10 B) \sec ^3(c+d x) \tan (c+d x)}{63 d \sqrt {a+a \sec (c+d x)}}+\frac {2 a B \sec ^3(c+d x) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{9 d}+\frac {1}{21} (a (39 A+34 B)) \int \sec ^3(c+d x) \sqrt {a+a \sec (c+d x)} \, dx \\ & = \frac {2 a^2 (9 A+10 B) \sec ^3(c+d x) \tan (c+d x)}{63 d \sqrt {a+a \sec (c+d x)}}+\frac {2 a B \sec ^3(c+d x) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{9 d}+\frac {2 (39 A+34 B) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}+\frac {1}{105} (2 (39 A+34 B)) \int \sec (c+d x) \left (\frac {3 a}{2}-a \sec (c+d x)\right ) \sqrt {a+a \sec (c+d x)} \, dx \\ & = \frac {2 a^2 (9 A+10 B) \sec ^3(c+d x) \tan (c+d x)}{63 d \sqrt {a+a \sec (c+d x)}}-\frac {4 a (39 A+34 B) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a B \sec ^3(c+d x) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{9 d}+\frac {2 (39 A+34 B) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d}+\frac {1}{45} (a (39 A+34 B)) \int \sec (c+d x) \sqrt {a+a \sec (c+d x)} \, dx \\ & = \frac {2 a^2 (39 A+34 B) \tan (c+d x)}{45 d \sqrt {a+a \sec (c+d x)}}+\frac {2 a^2 (9 A+10 B) \sec ^3(c+d x) \tan (c+d x)}{63 d \sqrt {a+a \sec (c+d x)}}-\frac {4 a (39 A+34 B) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{315 d}+\frac {2 a B \sec ^3(c+d x) \sqrt {a+a \sec (c+d x)} \tan (c+d x)}{9 d}+\frac {2 (39 A+34 B) (a+a \sec (c+d x))^{3/2} \tan (c+d x)}{105 d} \\
\end{align*}
Mathematica [A] (verified)
Time = 0.49 (sec) , antiderivative size = 100, normalized size of antiderivative = 0.53
\[
\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx=\frac {2 a^2 \left (8 (39 A+34 B)+4 (39 A+34 B) \sec (c+d x)+3 (39 A+34 B) \sec ^2(c+d x)+5 (9 A+17 B) \sec ^3(c+d x)+35 B \sec ^4(c+d x)\right ) \tan (c+d x)}{315 d \sqrt {a (1+\sec (c+d x))}}
\]
[In]
Integrate[Sec[c + d*x]^3*(a + a*Sec[c + d*x])^(3/2)*(A + B*Sec[c + d*x]),x]
[Out]
(2*a^2*(8*(39*A + 34*B) + 4*(39*A + 34*B)*Sec[c + d*x] + 3*(39*A + 34*B)*Sec[c + d*x]^2 + 5*(9*A + 17*B)*Sec[c
+ d*x]^3 + 35*B*Sec[c + d*x]^4)*Tan[c + d*x])/(315*d*Sqrt[a*(1 + Sec[c + d*x])])
Maple [A] (verified)
Time = 4.30 (sec) , antiderivative size = 131, normalized size of antiderivative = 0.69
| | |
| method | result | size |
| | |
| default |
\(\frac {2 a \left (312 A \cos \left (d x +c \right )^{4}+272 B \cos \left (d x +c \right )^{4}+156 A \cos \left (d x +c \right )^{3}+136 B \cos \left (d x +c \right )^{3}+117 A \cos \left (d x +c \right )^{2}+102 B \cos \left (d x +c \right )^{2}+45 A \cos \left (d x +c \right )+85 B \cos \left (d x +c \right )+35 B \right ) \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \tan \left (d x +c \right ) \sec \left (d x +c \right )^{3}}{315 d \left (\cos \left (d x +c \right )+1\right )}\) |
\(131\) |
| parts |
\(\frac {2 A a \left (104 \cos \left (d x +c \right )^{3}+52 \cos \left (d x +c \right )^{2}+39 \cos \left (d x +c \right )+15\right ) \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \tan \left (d x +c \right ) \sec \left (d x +c \right )^{2}}{105 d \left (\cos \left (d x +c \right )+1\right )}+\frac {2 B a \left (272 \cos \left (d x +c \right )^{4}+136 \cos \left (d x +c \right )^{3}+102 \cos \left (d x +c \right )^{2}+85 \cos \left (d x +c \right )+35\right ) \sqrt {a \left (1+\sec \left (d x +c \right )\right )}\, \tan \left (d x +c \right ) \sec \left (d x +c \right )^{3}}{315 d \left (\cos \left (d x +c \right )+1\right )}\) |
\(158\) |
| | |
|
|
|
|
[In]
int(sec(d*x+c)^3*(a+a*sec(d*x+c))^(3/2)*(A+B*sec(d*x+c)),x,method=_RETURNVERBOSE)
[Out]
2/315*a/d*(312*A*cos(d*x+c)^4+272*B*cos(d*x+c)^4+156*A*cos(d*x+c)^3+136*B*cos(d*x+c)^3+117*A*cos(d*x+c)^2+102*
B*cos(d*x+c)^2+45*A*cos(d*x+c)+85*B*cos(d*x+c)+35*B)*(a*(1+sec(d*x+c)))^(1/2)/(cos(d*x+c)+1)*tan(d*x+c)*sec(d*
x+c)^3
Fricas [A] (verification not implemented)
none
Time = 0.25 (sec) , antiderivative size = 127, normalized size of antiderivative = 0.67
\[
\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx=\frac {2 \, {\left (8 \, {\left (39 \, A + 34 \, B\right )} a \cos \left (d x + c\right )^{4} + 4 \, {\left (39 \, A + 34 \, B\right )} a \cos \left (d x + c\right )^{3} + 3 \, {\left (39 \, A + 34 \, B\right )} a \cos \left (d x + c\right )^{2} + 5 \, {\left (9 \, A + 17 \, B\right )} a \cos \left (d x + c\right ) + 35 \, B a\right )} \sqrt {\frac {a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sin \left (d x + c\right )}{315 \, {\left (d \cos \left (d x + c\right )^{5} + d \cos \left (d x + c\right )^{4}\right )}}
\]
[In]
integrate(sec(d*x+c)^3*(a+a*sec(d*x+c))^(3/2)*(A+B*sec(d*x+c)),x, algorithm="fricas")
[Out]
2/315*(8*(39*A + 34*B)*a*cos(d*x + c)^4 + 4*(39*A + 34*B)*a*cos(d*x + c)^3 + 3*(39*A + 34*B)*a*cos(d*x + c)^2
+ 5*(9*A + 17*B)*a*cos(d*x + c) + 35*B*a)*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sin(d*x + c)/(d*cos(d*x + c)
^5 + d*cos(d*x + c)^4)
Sympy [F]
\[
\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx=\int \left (a \left (\sec {\left (c + d x \right )} + 1\right )\right )^{\frac {3}{2}} \left (A + B \sec {\left (c + d x \right )}\right ) \sec ^{3}{\left (c + d x \right )}\, dx
\]
[In]
integrate(sec(d*x+c)**3*(a+a*sec(d*x+c))**(3/2)*(A+B*sec(d*x+c)),x)
[Out]
Integral((a*(sec(c + d*x) + 1))**(3/2)*(A + B*sec(c + d*x))*sec(c + d*x)**3, x)
Maxima [F]
\[
\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx=\int { {\left (B \sec \left (d x + c\right ) + A\right )} {\left (a \sec \left (d x + c\right ) + a\right )}^{\frac {3}{2}} \sec \left (d x + c\right )^{3} \,d x }
\]
[In]
integrate(sec(d*x+c)^3*(a+a*sec(d*x+c))^(3/2)*(A+B*sec(d*x+c)),x, algorithm="maxima")
[Out]
8/315*(315*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)*(((3*A + 2*B)*a*d*cos(2*d*
x + 2*c)^4 + (3*A + 2*B)*a*d*sin(2*d*x + 2*c)^4 + 4*(3*A + 2*B)*a*d*cos(2*d*x + 2*c)^3 + 6*(3*A + 2*B)*a*d*cos
(2*d*x + 2*c)^2 + 4*(3*A + 2*B)*a*d*cos(2*d*x + 2*c) + (3*A + 2*B)*a*d + 2*((3*A + 2*B)*a*d*cos(2*d*x + 2*c)^2
+ 2*(3*A + 2*B)*a*d*cos(2*d*x + 2*c) + (3*A + 2*B)*a*d)*sin(2*d*x + 2*c)^2)*integrate((cos(2*d*x + 2*c)^2 + s
in(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)*(((cos(10*d*x + 10*c)*cos(2*d*x + 2*c) + 4*cos(8*d*x + 8*c)*
cos(2*d*x + 2*c) + 6*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 4*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c
)^2 + sin(10*d*x + 10*c)*sin(2*d*x + 2*c) + 4*sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 6*sin(6*d*x + 6*c)*sin(2*d*x
+ 2*c) + 4*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(7/2*arctan2(sin(2*d*x + 2*c), cos(2*d*
x + 2*c))) + (cos(2*d*x + 2*c)*sin(10*d*x + 10*c) + 4*cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 6*cos(2*d*x + 2*c)*s
in(6*d*x + 6*c) + 4*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(10*d*x + 10*c)*sin(2*d*x + 2*c) - 4*cos(8*d*x + 8*
c)*sin(2*d*x + 2*c) - 6*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 4*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(7/2*arcta
n2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*
x + 2*c)*sin(10*d*x + 10*c) + 4*cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 6*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 4*co
s(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(10*d*x + 10*c)*sin(2*d*x + 2*c) - 4*cos(8*d*x + 8*c)*sin(2*d*x + 2*c) -
6*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 4*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(7/2*arctan2(sin(2*d*x + 2*c), c
os(2*d*x + 2*c))) - (cos(10*d*x + 10*c)*cos(2*d*x + 2*c) + 4*cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 6*cos(6*d*x +
6*c)*cos(2*d*x + 2*c) + 4*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(10*d*x + 10*c)*sin(2*d
*x + 2*c) + 4*sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 6*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 4*sin(4*d*x + 4*c)*sin
(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*sin(7/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(3/2*arctan2(sin(
2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))/((cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*
d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*
cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)*cos(6*d*x + 6*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)
^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(10*d*x + 10
*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 36*(cos(2*d
*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2
*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*s
in(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(
2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos
(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos
(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(10*d*x + 10*c) + 8*(cos(2*d*x + 2*c)^3 + cos(2*d*
x + 2*c)*sin(2*d*x + 2*c)^2 + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x +
6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x +
2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 12*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4
*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 +
cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 8*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x
+ 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 4*(cos(2*d*x + 2
*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*
c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2
*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(10*d*x + 10*c)
+ 8*(sin(2*d*x + 2*c)^3 + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*
c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)
^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 12*(sin(2*d*x + 2*c)^3 + 4*(cos(2*d*x + 2*c)
^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 8*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1
)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x
+ 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(10*d*x
+ 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 36*(cos
(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 16*(cos(2*d*x + 2*c)^2 + s
in(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 +
sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^
2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2
*c) + 1)*sin(6*d*x + 6*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x
+ 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*
d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x
+ 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x +
2*c))*cos(10*d*x + 10*c) + 8*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 6*(cos(2*d*x + 2*c)^
2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2
+ 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 12*(
cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(
2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 8*(cos(2*d*x
+ 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + c
os(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)*sin(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) +
4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(10*d*x + 10*c) + 8*(sin(2*d*x + 2*c)^3 + 6*(cos(2*d*x + 2*c)^2 +
sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin
(8*d*x + 8*c) + 12*(sin(2*d*x + 2*c)^3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*
sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 8*(sin(2
*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(3/2*arct
an2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2), x) + 2*((2*A + 3*B)*a*d*cos(2*d*x + 2*c)^4 + (2*A + 3*B)*a*d*
sin(2*d*x + 2*c)^4 + 4*(2*A + 3*B)*a*d*cos(2*d*x + 2*c)^3 + 6*(2*A + 3*B)*a*d*cos(2*d*x + 2*c)^2 + 4*(2*A + 3*
B)*a*d*cos(2*d*x + 2*c) + (2*A + 3*B)*a*d + 2*((2*A + 3*B)*a*d*cos(2*d*x + 2*c)^2 + 2*(2*A + 3*B)*a*d*cos(2*d*
x + 2*c) + (2*A + 3*B)*a*d)*sin(2*d*x + 2*c)^2)*integrate((cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d
*x + 2*c) + 1)^(1/4)*(((cos(10*d*x + 10*c)*cos(2*d*x + 2*c) + 4*cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 6*cos(6*d*
x + 6*c)*cos(2*d*x + 2*c) + 4*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(10*d*x + 10*c)*sin(
2*d*x + 2*c) + 4*sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 6*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 4*sin(4*d*x + 4*c)*
sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) + (cos(2*d*x + 2*c
)*sin(10*d*x + 10*c) + 4*cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 6*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 4*cos(2*d*x
+ 2*c)*sin(4*d*x + 4*c) - cos(10*d*x + 10*c)*sin(2*d*x + 2*c) - 4*cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 6*cos(6
*d*x + 6*c)*sin(2*d*x + 2*c) - 4*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*
x + 2*c))))*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(10*d*x + 10*c) +
4*cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 6*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 4*cos(2*d*x + 2*c)*sin(4*d*x + 4*
c) - cos(10*d*x + 10*c)*sin(2*d*x + 2*c) - 4*cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 6*cos(6*d*x + 6*c)*sin(2*d*x
+ 2*c) - 4*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) - (cos(10*d
*x + 10*c)*cos(2*d*x + 2*c) + 4*cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 6*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 4*co
s(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(10*d*x + 10*c)*sin(2*d*x + 2*c) + 4*sin(8*d*x + 8*c
)*sin(2*d*x + 2*c) + 6*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 4*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2
*c)^2)*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c
) + 1)))/((cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x +
2*c) + 1)*cos(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d
*x + 8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 16*(c
os(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (
cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^
2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c
)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x +
2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*
d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x
+ 2*c) + 1)*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x
+ 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x +
2*c)^2 + cos(2*d*x + 2*c))*cos(10*d*x + 10*c) + 8*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 +
6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2
+ sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*co
s(8*d*x + 8*c) + 12*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(2*
d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x +
6*c) + 8*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))
*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1
)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (
cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(10*d*x + 10*c) + 8*(sin(2*d*x + 2*c)^3 + 6*
(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 +
sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*
sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 12*(sin(2*d*x + 2*c)^3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*c
os(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*
d*x + 6*c) + 8*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x
+ 4*c))*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4
+ (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2
*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d
*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*
d*x + 2*c) + 1)*sin(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*s
in(8*d*x + 8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 +
16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c
)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2
+ 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)
^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^
2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(10*d*x + 10*c) + 8
*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*co
s(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*co
s(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 12*(cos(2*d*x + 2*c)^3 + cos(2*d*
x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x +
4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 8*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*s
in(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*
d*x + 2*c)^3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c) + 6*(cos(
2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2
*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2
*d*x + 2*c))*sin(10*d*x + 10*c) + 8*(sin(2*d*x + 2*c)^3 + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2
*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4
*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 12*(sin(2*d*x
+ 2*c)^3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x
+ 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 8*(sin(2*d*x + 2*c)^3 + (cos(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d
*x + 2*c) + 1))^2), x) + (A*a*d*cos(2*d*x + 2*c)^4 + A*a*d*sin(2*d*x + 2*c)^4 + 4*A*a*d*cos(2*d*x + 2*c)^3 + 6
*A*a*d*cos(2*d*x + 2*c)^2 + 4*A*a*d*cos(2*d*x + 2*c) + A*a*d + 2*(A*a*d*cos(2*d*x + 2*c)^2 + 2*A*a*d*cos(2*d*x
+ 2*c) + A*a*d)*sin(2*d*x + 2*c)^2)*integrate((cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)^(1/4)*(((cos(10*d*x + 10*c)*cos(2*d*x + 2*c) + 4*cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 6*cos(6*d*x + 6*c)*co
s(2*d*x + 2*c) + 4*cos(4*d*x + 4*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(10*d*x + 10*c)*sin(2*d*x + 2*c
) + 4*sin(8*d*x + 8*c)*sin(2*d*x + 2*c) + 6*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 4*sin(4*d*x + 4*c)*sin(2*d*x +
2*c) + sin(2*d*x + 2*c)^2)*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) + (cos(2*d*x + 2*c)*sin(10*d*
x + 10*c) + 4*cos(2*d*x + 2*c)*sin(8*d*x + 8*c) + 6*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 4*cos(2*d*x + 2*c)*sin
(4*d*x + 4*c) - cos(10*d*x + 10*c)*sin(2*d*x + 2*c) - 4*cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 6*cos(6*d*x + 6*c)
*sin(2*d*x + 2*c) - 4*cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))
*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((cos(2*d*x + 2*c)*sin(10*d*x + 10*c) + 4*cos(2*d*
x + 2*c)*sin(8*d*x + 8*c) + 6*cos(2*d*x + 2*c)*sin(6*d*x + 6*c) + 4*cos(2*d*x + 2*c)*sin(4*d*x + 4*c) - cos(10
*d*x + 10*c)*sin(2*d*x + 2*c) - 4*cos(8*d*x + 8*c)*sin(2*d*x + 2*c) - 6*cos(6*d*x + 6*c)*sin(2*d*x + 2*c) - 4*
cos(4*d*x + 4*c)*sin(2*d*x + 2*c))*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))) - (cos(10*d*x + 10*c)*
cos(2*d*x + 2*c) + 4*cos(8*d*x + 8*c)*cos(2*d*x + 2*c) + 6*cos(6*d*x + 6*c)*cos(2*d*x + 2*c) + 4*cos(4*d*x + 4
*c)*cos(2*d*x + 2*c) + cos(2*d*x + 2*c)^2 + sin(10*d*x + 10*c)*sin(2*d*x + 2*c) + 4*sin(8*d*x + 8*c)*sin(2*d*x
+ 2*c) + 6*sin(6*d*x + 6*c)*sin(2*d*x + 2*c) + 4*sin(4*d*x + 4*c)*sin(2*d*x + 2*c) + sin(2*d*x + 2*c)^2)*sin(
3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c))))*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))/((
cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*c
os(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2
+ 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 16*(cos(2*d*x +
2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x +
2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d
*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos
(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*
sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^
3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)
*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*
(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + c
os(2*d*x + 2*c))*cos(10*d*x + 10*c) + 8*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 6*(cos(2*d
*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*
x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8
*c) + 12*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^
2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 8*(
cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x
+ 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x
+ 2*c) + 1)*sin(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x
+ 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(10*d*x + 10*c) + 8*(sin(2*d*x + 2*c)^3 + 6*(cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x +
2*c))*sin(8*d*x + 8*c) + 12*(sin(2*d*x + 2*c)^3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x +
2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c)
+ 8*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*co
s(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))^2 + (cos(2*d*x + 2*c)^4 + sin(2*d*x + 2*c)^4 + (cos(2*d
*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin
(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2
*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) +
1)*cos(4*d*x + 4*c)^2 + 2*cos(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(10*d*x + 10*c)^2 + 16*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x +
8*c)^2 + 36*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c)^2 + 16*(cos(2*
d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)^2 + (2*cos(2*d*x + 2*c)^2 + 2*cos
(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)^2 + 2*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 4*(cos(2
*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c)^2 + sin(2*
d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2
*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(10*d*x + 10*c) + 8*(cos(2*d*x
+ 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x + 2*c)^2 + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2
*c) + 1)*cos(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4
*c) + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(8*d*x + 8*c) + 12*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*si
n(2*d*x + 2*c)^2 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*cos(4*d*x + 4*c) + 2*c
os(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(6*d*x + 6*c) + 8*(cos(2*d*x + 2*c)^3 + cos(2*d*x + 2*c)*sin(2*d*x +
2*c)^2 + 2*cos(2*d*x + 2*c)^2 + cos(2*d*x + 2*c))*cos(4*d*x + 4*c) + cos(2*d*x + 2*c)^2 + 2*(sin(2*d*x + 2*c)^
3 + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(8*d*x + 8*c) + 6*(cos(2*d*x + 2*c
)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)
^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c)
)*sin(10*d*x + 10*c) + 8*(sin(2*d*x + 2*c)^3 + 6*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c)
+ 1)*sin(6*d*x + 6*c) + 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c)
+ (cos(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(8*d*x + 8*c) + 12*(sin(2*d*x + 2*c)^3 +
4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)*sin(4*d*x + 4*c) + (cos(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(6*d*x + 6*c) + 8*(sin(2*d*x + 2*c)^3 + (cos(2*d*x + 2*c)^2 + 2*
cos(2*d*x + 2*c) + 1)*sin(2*d*x + 2*c))*sin(4*d*x + 4*c))*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) +
1))^2), x))*sqrt(a) - (3*(35*A*a*sin(6*d*x + 6*c) + 126*(A + B)*a*sin(4*d*x + 4*c) + 3*(39*A + 34*B)*a*sin(2*
d*x + 2*c))*cos(9/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - (105*A*a*cos(6*d*x + 6*c) + 378*(A + B)
*a*cos(4*d*x + 4*c) + 9*(39*A + 34*B)*a*cos(2*d*x + 2*c) + 2*(39*A + 34*B)*a)*sin(9/2*arctan2(sin(2*d*x + 2*c)
, cos(2*d*x + 2*c) + 1)))*sqrt(a))/((d*cos(2*d*x + 2*c)^4 + d*sin(2*d*x + 2*c)^4 + 4*d*cos(2*d*x + 2*c)^3 + 6*
d*cos(2*d*x + 2*c)^2 + 2*(d*cos(2*d*x + 2*c)^2 + 2*d*cos(2*d*x + 2*c) + d)*sin(2*d*x + 2*c)^2 + 4*d*cos(2*d*x
+ 2*c) + d)*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4))
Giac [F]
\[
\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx=\int { {\left (B \sec \left (d x + c\right ) + A\right )} {\left (a \sec \left (d x + c\right ) + a\right )}^{\frac {3}{2}} \sec \left (d x + c\right )^{3} \,d x }
\]
[In]
integrate(sec(d*x+c)^3*(a+a*sec(d*x+c))^(3/2)*(A+B*sec(d*x+c)),x, algorithm="giac")
[Out]
sage0*x
Mupad [B] (verification not implemented)
Time = 22.81 (sec) , antiderivative size = 596, normalized size of antiderivative = 3.15
\[
\int \sec ^3(c+d x) (a+a \sec (c+d x))^{3/2} (A+B \sec (c+d x)) \, dx=\frac {\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (-\frac {a\,\left (3\,A+2\,B\right )\,8{}\mathrm {i}}{7\,d}+\frac {a\,\left (A+4\,B\right )\,8{}\mathrm {i}}{7\,d}+\frac {B\,a\,32{}\mathrm {i}}{63\,d}\right )+\frac {A\,a\,8{}\mathrm {i}}{7\,d}-\frac {a\,\left (A+2\,B\right )\,24{}\mathrm {i}}{7\,d}-\frac {B\,a\,32{}\mathrm {i}}{7\,d}\right )}{\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )\,{\left ({\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}+1\right )}^3}+\frac {\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (-\frac {a\,\left (2\,A+3\,B\right )\,16{}\mathrm {i}}{9\,d}+\frac {a\,\left (3\,A+2\,B\right )\,8{}\mathrm {i}}{9\,d}+\frac {A\,a\,8{}\mathrm {i}}{9\,d}\right )+\frac {a\,\left (2\,A+3\,B\right )\,16{}\mathrm {i}}{9\,d}-\frac {a\,\left (3\,A+2\,B\right )\,8{}\mathrm {i}}{9\,d}-\frac {A\,a\,8{}\mathrm {i}}{9\,d}\right )}{\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )\,{\left ({\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}+1\right )}^4}+\frac {\left (\frac {A\,a\,8{}\mathrm {i}}{3\,d}-\frac {a\,{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (39\,A+34\,B\right )\,8{}\mathrm {i}}{315\,d}\right )\,\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}}{\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )\,\left ({\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}+1\right )}+\frac {\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\left (\frac {a\,\left (3\,A+2\,B\right )\,8{}\mathrm {i}}{5\,d}+\frac {a\,\left (3\,A+B\right )\,16{}\mathrm {i}}{105\,d}\right )-\frac {A\,a\,8{}\mathrm {i}}{5\,d}+\frac {a\,\left (A+3\,B\right )\,16{}\mathrm {i}}{5\,d}\right )}{\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )\,{\left ({\mathrm {e}}^{c\,2{}\mathrm {i}+d\,x\,2{}\mathrm {i}}+1\right )}^2}-\frac {a\,{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}\,\sqrt {a+\frac {a}{\frac {{\mathrm {e}}^{-c\,1{}\mathrm {i}-d\,x\,1{}\mathrm {i}}}{2}+\frac {{\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}}{2}}}\,\left (39\,A+34\,B\right )\,16{}\mathrm {i}}{315\,d\,\left ({\mathrm {e}}^{c\,1{}\mathrm {i}+d\,x\,1{}\mathrm {i}}+1\right )}
\]
[In]
int(((A + B/cos(c + d*x))*(a + a/cos(c + d*x))^(3/2))/cos(c + d*x)^3,x)
[Out]
((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2)*(exp(c*1i + d*x*1i)*((a*(A + 4*B)*8i)/(7*d) - (
a*(3*A + 2*B)*8i)/(7*d) + (B*a*32i)/(63*d)) + (A*a*8i)/(7*d) - (a*(A + 2*B)*24i)/(7*d) - (B*a*32i)/(7*d)))/((e
xp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)^3) + ((a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^
(1/2)*(exp(c*1i + d*x*1i)*((a*(3*A + 2*B)*8i)/(9*d) - (a*(2*A + 3*B)*16i)/(9*d) + (A*a*8i)/(9*d)) + (a*(2*A +
3*B)*16i)/(9*d) - (a*(3*A + 2*B)*8i)/(9*d) - (A*a*8i)/(9*d)))/((exp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) +
1)^4) + (((A*a*8i)/(3*d) - (a*exp(c*1i + d*x*1i)*(39*A + 34*B)*8i)/(315*d))*(a + a/(exp(- c*1i - d*x*1i)/2 + e
xp(c*1i + d*x*1i)/2))^(1/2))/((exp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)) + ((a + a/(exp(- c*1i - d*x*1
i)/2 + exp(c*1i + d*x*1i)/2))^(1/2)*(exp(c*1i + d*x*1i)*((a*(3*A + 2*B)*8i)/(5*d) + (a*(3*A + B)*16i)/(105*d))
- (A*a*8i)/(5*d) + (a*(A + 3*B)*16i)/(5*d)))/((exp(c*1i + d*x*1i) + 1)*(exp(c*2i + d*x*2i) + 1)^2) - (a*exp(c
*1i + d*x*1i)*(a + a/(exp(- c*1i - d*x*1i)/2 + exp(c*1i + d*x*1i)/2))^(1/2)*(39*A + 34*B)*16i)/(315*d*(exp(c*1
i + d*x*1i) + 1))